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ABSTRACT 

In this paper, we model dwarf galaxies as a two-component system of grav- 
itationally coupled stars and atomic hydrogen gas in the external force field 
of a pseudo-isothermal dark matter halo, and numerically obtain the radial 
distribution of H i vertical scale heights. This is done for a group of four dwarf 
galaxies (DDO 154, Ho II, IC 2574 and NGC 2366) for which most necessary 
input parameters are available from observations. The formulation of the equa- 
tions takes into account the rising rotation curves generally observed in dwarf 
galaxies. The inclusion of self-gravity of the gas into the model at par with 
that of the stars results in scale heights that are smaller than what was ob- 
tained by previous authors. This is important as the gas scale height is often 
used for deriving other physical quantities. The inclusion of gas self-gravity 
is particularly relevant in the case of dwarf galaxies where the gas cannot be 
considered a minor perturbation to the mass distribution of the stars. We find 
that three out of four galaxies studied show a flaring of their Hi disks with 
increasing radius, by a factor of a few within several disk scale lengths. The 
fourth galaxy has a thick Hi disk throughout. This arises as a result of the 
gas velocity dispersion remaining constant or decreasing only slightly while 
the disk mass distribution declines exponentially as a function of radius. 
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1 INTRODUCTION 



In this paper we address the vertical distribution of neutral atomic hydrogen (Hi) in four 
ar (dl rr) galaxies that are part of THINGS (The Hi Nearby Galaxy Survey; 



dwarf irregu 


ar (dl 


Walter et al. 


2008) 



20081 ). These objects occupy the low mass, low luminosity tail of the galaxy 



luminosity function. They have accordingly low metallicities, in agreement with the find- 
ing that star formation in these systems occurs episodically, in bursts, with long periods 
during which apparently little star formation is taking place. Contrary to larger spiral 
galaxies, where stars dominate, in dwarf galaxies stars and Hi make up the baryonic mass 
content on an equal footing, with little or no H2 being trac ed through CO observations 



( iTacconi fc Young] Il987l : iTaylor. Kobulnicky. fc Skillmaru 119981 ). The jury is still out on the 



relation of CO to H2 in low-metallicity, low-mass environments encountered in dwarf galax- 



ies, although it seems clear that a low CO signa 



low H2 mass (jWilson 



1995 



Bolatto et al. 



2008 



does not necessar ily mean a correspondingly 



Leroy et al.ll2009l ). Their low metallicity and 



the fact that a large fraction of their gas has not yet been turned into stars is often used to ar- 
gue that dwarf galaxies, within the paradigm of hierarchical galaxy formation, then must be 
broadly similar to the building blocks at high redshift that eventually formed today's larger 
galaxies. In accordance with models of hierarchical structure formati on dwarf galaxies are 



dark matter domina t ed, something w hich is born out by observations (ICarignan &: Freeman 



1988 



Ashman 



1992 



Oh et al.l 



20101). They co nstitute the largest number of galaxies in the 



present day observable universe ( lMateolll998l ). 



Brinks. Walter. &: OttI (|2002[ ) postulated that the gas scale height in dirrs was larger, 
both in absolute as well as relative sense, than in larger, spiral galaxies. This has several 
implications, especially if indeed nearby dIrrs are broadly similar to those found at large 
look-back times. For example, their cross section for intercepting distant QSOs, giving rise 
to damped Lyman-a signatures, will be much larger. Also, a thicker gas layer will tend to 
allow Hi shells and bubbles to grow to larger diameters, preventing early break-out and 
thus reducing the likelihood that metals pollute the intergalactic medium. 

The vertical stellar distri bution has been ext ensively studied observationally (see for 



a review of his early papers 



van der Kruit 



19881 ). The fact that t he stellar distr i butio n 



shows a moderate flaring at large radii was shown obser vationallv bv ide 



and predicted from a general self-consistent model by 



Griis et all (Il997f ). 



Naravan k. Jog! (l2002bl ). In broad 



terms, the observed constant velocity dispersion of the gas combined with a decreasing 
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mass surface density (as traced by the exponential decline of the stellar light distribution) 
predi cts a flaring Hi layer. Early direct observation s of this flaring were m ade in the Milky 



Way (iHenderson. Jackson, fc Kerr 



1982) and M31 flBrinks fc Burtonlll984h 



Observing the Hi gas scale height directly is difficult, even for a perfectly edge-on 
galaxy, due to effects such as line-of -sight integration , beam smearing and optical depth 
effects as was already point ed out bv ISancisi fc AllenI (119791). More recent observation s of 



the Milky Way are given by 



also 



Wouterloot et a. 



Kalberla fc KerpI 120091 ). and for M31 by 



199o[ ) 



Braun 



and 



Levine. Blitz, fc HeilesI ( 12006 



see 



(Il99ll ). S ome other syst ems whose scale 



heights have been measured are NGC 3198 and NGC 2403 (jSickine] 119971 ). No serious at- 
tempt has been made to date to investigate th e vertical stellar and gas distributions in dirrs. 

(l2002[ l. because their observed Hi velocity dispersion 



As discussed by 



Brinks. Walter, fc Ott 



is not very different from that observed in spirals, but their mass surface density is consid- 
erably lower, dIrrs must have a larger scale height in Hi as compared to spiral galaxies in 
an absolute sense. 

The origin of flaring of the gas layer in spiral g alaxies has been treated recently in a 
self-consistently manner by iNarayan &: Jog) (l2002al l. who produced a numerical model of 
the vertical scale heights of the gas and stars in a gravitationally-coupled, three-component 
Galactic disk. In this paper, we extend this earlier work on the vertical stellar and gas 
distribution in spirals to dwarf galaxies, and determine the scale height of the g 
function of radius analytically. We model the galaxy as a gravitationally coupled system of 
stars and gas subjected to the external force field of a dark matter halo, and use the model 
to obtain the H I vertical scale height as a function of galacto-centr ic radius for four dwarf 



galaxies. In earlier papers we had assumed a fiat rotation curve (INarayan &: Jog 



Narayan. Saha. fc Jog 



2002 



3 



20051 ). This approximation is not warranted for dwarf galaxies which 



show a steadily rising rotation curve throughout most of their disk. Instead we explicitly 
take account here of a rising rotation curve, and show that the originally two-dimensional 
partial differential equations can still be written as one-dimensional or ordinary differential 
equations, as in the case for the giant spiral galaxies which have flat rotation curves. 

We will concentrate in this paper on the properties of the bulk of the gas and therefore 
ignore the con tribution by the ~ 10 percent of "extra-planar gas" fo und by Fraternali and 



collaborators ( iMarasco &: FraternaU 



2011 



Fraternali Sz Binney 



20081 ) . most of which seems 



to be related to energy input from supernovae into the interstellar medium. We describe our 



4 A. Banerjee et al. 

methodology in Sects. [2] and [3l justify the values for the input parameters to the model in 
Sect, m and present our results and their discussion in Sects. [5l 



2 MODEL 



2.1 Gravitationally— coupled, two— component, galactic disk model 

We use a two-component galactic disk model of gravitationally-coupled stars and Hi gas 
subjected to the external force-field of a dark matter halo with a pseudo-isothermal den- 
sity profile. In this model, the stars and gas are assumed to be present in the form of 



concentric, thin, axis ymmetric disks embedded within each other ( Narayan &: Jo 



Banerjee fc JoglbnOTl ). 



2002b : 



The joint Poisson equation for an axisymmetric system of stars and gas is given by: 



total ' 



RdR^ OR 



+ 



ATiG{Y,Pi + Ph) 



(1) 



i=l 



where iptotai is the net potential of the disk due to the stars, the H I gas and the dark matter 
halo, pi with z = 1 to 2 denotes the mass volume density for each of the disk components 
(star s, gas) while ph denotes the same for the halo, and is given by (Ide Zeeuw fc Pfenniger 
lOSsI): 



ph{R,z) 



Po 



1 



(2) 



R'i 



where po is the central core density of the halo, and Rc is the core radius, as apphcable 
for a pseudo-isothermal density distribution. One may note here that since dwarf galaxies 
have rotation curves that are generally rising out to the last measured point, the radial term 
in equa tion (1) cannot be ignored as in the case for large spirals with almost flat rotation 



curves (INarayan. Saha. fc Jog 



2005 



Baneriee fc Jogll2008h . 



The e quatiori of hy drostatic equilibrium for the i^^ disk component in the z direction is 
given by (IRohlfsHlOTTf ): 



total 







(3) 



where < v"^ >i is the mean square velocity of the i*^^ disk component, equation (3) 
constitutes two equations, one each for the density of each of the two disk components. In 
order to obtain the vertical density distribution and hence the scale height of Hi at each 
R, one needs to solve equation (1) and the set of two equations given by equation 3 jointly, 
which together constitute a set of three coupled differential equations in the variables pi and 
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4'totai in two dimensions R, z. Even if all the initial conditions are known, obtaining a rigorous 
solution of the above equations is going to be extremely complicated and computationally 
expensive. Thus we look for suitable approximations to simplify the problem. 

2.2 Approximations 

Combining equation (1) and (3) and assuming isothermal condition for each of the two disk 
components we get: 

Now, we replace the radial derivative term by its approximate value by utilising our 
knowledge of the shape of the observed rotation curve, which will effectively reduce the 
above equation to a one dimensional problem along the z direction, as follows: 

where {vrot)Rz intrinsic rotational velocity at any R,z. However, we do not a 

priori know {vrot)R,z as a function of 2; at a given R. However, at each i?, the observed 
rotation curve gives the intensity-weighted average of {vrot)R,z along the vertical direction. 
We approximate the value of (frot)_R,z at each z by the commonly used intensity-weighted 
average velocity, i.e., the rotational velocity given by the observed rotation curve. As a result 
of this approximation, equation (4) reduces to 

The last term on the r.h.s now becomes a function of R only, and can be easily determined 
from the gradient of the observed rotation curve. We note that the presence of the last 
term on the r.h.s dilutes the effect of the net local mass density as far as the downward 
gravitational pull is concerned. Simple calculations show that in the rising part of the rotation 
curve the value of the term is of the same order of magnitude, but entering with opposite 
sign, as that of the net local gravity which results in a more puffed-up gas distribution and 
hence a larger value for the gas scale height. 

We further stress that in solving equation (6), we have used the dark matter core density 
(p/i) and the last term involving the gradient of the rotational velocity, both of which have 
been obtained from the observed rotation curve. While the latter has been determined 
directly from the observed curve, was estimated after subtracting the disk contribution 
to the rotational velocity based on the calculated mass model of the disk. 
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3 NUMERICAL CALCULATIONS 

The two coupled, second-order, ordinary differential equations given by equation (6) are 



solved nume rically using the 



fashion (see 



Narayan fc Jog 



burth order Runge-Kutta M ethod of integration in an iterative 



2002b 



Banerjee fc Jog 



20071 ). Since each of the two equations 



is a second-order differential equation, we need to state the following two initial conditions 
at the mid-plane (i.e., 2; = 0) for each component for solving the equations: 



dpi 

{po)i -1- 
az 







(7) 



However, the modified mid-plane density (po)i for each component is not known a priori. 
Instead the net surface density Sj(i?), given by twice the area under the curve of Pi{z) versus 
z, known observationally, is used as a secondary boundary condition. The required value of 
{pi)o is then determined by trial and error, and eventually the Pi{z) distribution is fixed. 



4 INPUT PARAMETERS 



As mentioned above, we have to provide as inputs the surface density S(i?) and the vertical 
velocity dispersion {a^)- =< {v'^)i >^^^ of each of the two components (stars and H l) in order 
to solve equation (6) at a given radius R, in addition to the observed r otation curve of the 



galaxy. Most of the input was taken from the compilation of THINGS 
papers. The stellar surface densities S^ (-R) are t a ken fr om 



except NGC2366 which is taken from 



Leroy et al. 



Walter et al 



20081) 



(120081 ) for all dwarfs 



Oh et all (I2OO8I ). Since a direct measure ment of the 



veloc ity dispersion of stars (cr*) is not available, it is calculated (as described in 



2008 



Appendix B) based on the scale length of the disk R^) and the surface density of 
the stars S*(i?). This assumes a stellar disk in vertical hydrostatic equilibrium, neglecting 
gas gravity and thus slightly underestimates the stellar velocity dispersion. However, the 
response of the gas to the net galactic potential is governed by the value of its own velocity 
dispersion (see equation (6)), and so a small change in the vertical density distribution of the 
stellar component can change the calc ulated gas s c ale he ights but marginally. The stellar 



disk scale lengths R^ are again from 



Leroy et al 



(I2OO8I ). We furthermore need the scale 



height of the stars (/i*) in order to calculate the volume density, (pn)*. W e used = Rd/7.3 
and an average flattening ratio (IKregel. van der Kruit. fc de Grijsll2002l ). 



We derived the Hi surfa ce density S^j f 



{(J^^) based on maps from 



Walter et al. 



R) a nd adopted a constant velocity dispersion 



(I2OO8I ). The resulting Hi scale height depends 
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strongly on the assumed value of the vertical velocity dispersion a^^ , which lies between 7 
- 9 kms~^. Our calculations show that a change in the value of the assumed Hi velocity 
dispersion by 1 kms^^ (or approximately 10 percent) changes the Hi scale height by 14-15 
percent. Therefore, the accuracy of the determined H I scale height depends on the accuracy 
with which a^^ can be determined. 



Tamburro et al. 



(120091), using the same THINGS data. 



published a more complete analysis of the velocity dispersions. Most galaxies show a gradual 
fall in Hi velocity dispersion with radius, ~ a few km s~^ per unit of R25. In addition to 
assuming a constant value, we use a smoothly falling velocity dispersion to model those two 
galaxies of our sample. Ho II and IC 2574, for which this more complete velocity dispersion 



information is available (ITamburro et al. 



20091 ). 



The dark matter central core densiti es Pn and core radii Rr are t aken from 



de Blok et al. 



( I2OO8I ) except for Ho II where we used [Bureau fc Carignaru (120021 ). These are determined 
from the observed rotation curves assuming rotational equilibrium, and a dark matter dis- 
tribution modelled a s a pseudo-isothermal halo. The ob served ro t ation curves are taken from 



de Blok et al. 



Bureau fc Carignan 



(12008 ) for D DQ 154 and NGC 2366, from lOh et all (12008[ l for IC 2574 and from 



( 120021 ) for Ho II. Key parameters used as input are hsted in Table 12] We 
note here that de Blok et al. (2008) estimated the disk contribution to the rotation curve by 
employing mass models for the stellar and the gaseous disks. This was subtracted from the 
observed rotation curve to obtain the contribution of the dark matter halo from which its 
core density and core radius are derived. The mass model of the disk is well constrained and 
allows for small variations in its parameters only that in turn may result in small changes 
in the deduced halo parameters; the effect on the gas scale height will be minimal. 

The essence of the radial mass distribution of the galaxy is effectively captured in the 
rotation curves. Dwarf galaxies in general have a rising rotation curve, i.e., the rotation 
velocity of these galaxies increases steadily with radius as far as 3-4 Rd, unlike spiral 
galaxies which sport flat rotation curves. We further note that in our previous work, we 
had modelled the observed gas scale height data to obt ain the best-fit halo parame ters, 



as was done for the Galaxv (iNaravan. Saha. fc Joe 



20051), M31 (IBaneriee fc JoeJ 120081 ) and 



UGC 7321 ( iBanerjee. Matthews. &: Jogll2010l ). Here the inverse approach is used, namely 



the observed rotation curve data is modelled to obtain the dark matter halo parameters and 
then the resulting H I vertical scale height distribution is derived. This approach is valid as 
long as the halo profile can be approximated by a pseudo-isothermal sphere. 
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5 RESULTS AND DISCUSSION 



We present in the top four panels of Figure [T] the resulting H I scale height versus radius for 
DDO 154, Ho II, IC2574 and NGC2366, assuming for each galaxy the constant Hi velocity 
dispersion as given in Tabled For each galaxy, the scale height calculations were done twice: 



first by neglecti ng the radial term in t 



spiral galaxies ( iNarayan fc Jog 



l e Poisson equation, m imicking the approach used for 



2nn2b 



Banerjee fc Jog) 120071 ). and a second time including 



the radial term. The two results differ from each other by about 10-20 percent, indicating 
that the exact shape of the rotation curve only enters as small, systematic correction. As 
mentioned above, the effect of the radial term is to somewhat counteract net local gravity, 
resulting in a larger value for the gas scale height. 

The individual data points follow, with some dispersion, a smooth increase with grow- 
ing radius. Most of this dispersion is due to radial variations in the input parameters and 
associated uncertainties, such as when calculating the derivative of the rotation curve. The 
combined effect of these uncertainties is of the order of 100 pc. 

A parameter which has a more substantial influence on the result is the velocity dis- 
persion, and its dependence on radius. The two bottom panels of Figure [1] show for two of 
the targets. Ho II and IC 2574, the effect of using a gradually declining velocity dispersion 
rather than a constant value. The rapid increase in H I scale height for IC 2574 becomes 
less pronounced as a result of the higher velocity dispersion near the centre, increasing the 
gas scale height there by about a factor of two, and a lower velocity dispersion near the 
last measured point, resulting in a relative increase across the gas disk of about a factor 
of two, rather than a factor of three earlier. The effect is even more pronounced for Ho II 
where we find an almost constant, but thick Hi disk throughout. This exercise illustrates 
the extent to which the resulting scale heights depend on accurate determinations of the 
velocity dispersion. 

In the remainder of this paper, we will use for all galaxies the analytical solution to the 
gas scale height while taking the radial term in the Poisson equation into account. In the case 
of Ho II and IC 2574 we take the the solution for the gradually decreasing velocity dispersion 
of the gas. Figure |2] then shows a pair of plots comparing the H I scale height distribution 
as a function of normalised radius. We use t wo normalisations, the first one by maximum 



radial extent of the Hi disc, -Rmax (taken from iBagetakos et al.ll2010l ). We also normalise by 



Gas scale heights in four dwarf galaxies 9 

Rsd, which corresponds to the radius where the mass surface density drops below 4 Mq pc~^. 
The -Rmax and the -Rsd values for the different galaxies have been listed in Table O 

We note that, with the exception of Ho II which has a thick gas disk throughout, the 
dwarfs all show H I flaring. The rate of flaring is highest in NGC 2366 and DDO 154 which 
can be ascribed to these two galaxies having a more centrally concentrated dark matter halo 
than the other two systems. 

Figure |3] presents a collection of plots of the logarithm of the scale height versus ra- 
dius, showing a near exponential behaviour, i.e. a linear slope, in each case. The curves for 
DDO 154 and IC 2574 display a break, necessitating two different slopes across two radial 
regimes, i.e., the gas scale height increases exponentially with radius with different scale 
lengths across the two regimes. It is interesting to compare this figure with Figure H] which 
shows the disk surface density of gas and stars combined versus radius. Knees in the graphs 
of DDO 154 and IC 2574 occur at the same radius as the breaks in the exponential curves 
in Fig. 131 confirming that the gas scale height depends strongly on the radial disk surface 
density distribution. 

This should make us pause for a moment as in the adopted radial surface density distri- 
bution distribution we have ignored any contribution by molecular gas. In the introduction 
we mentioned that dirr galaxies are deficient in CO which is used as a proxy to trace molec- 
ular gas. In fact, i n galaxies with metallicities below [12 + logO/H] < 7.5 no CO emission 
has been detected (ITaylor. Kobulnicky. fc Skillmarull998l ). We also mentioned in the intro- 
duction that a lack of CO does not necessarily imply a lac k of molecular gas. If present. 



this molecular gas 



2008 



Leroy et al. 



likely follows the stellar light distribution (lYoung et al. 



1995 



Bigiel et al. 



20091 ). Its effects on the flaring Hi beyond the stellar disks of the dwarfs 
will therefore be negligible. In the inner parts the molecular gas, if present, will contribute 
to the mass density, effectively reducing the scale height locally. 

Despite some modest potential adjustment at small radii due to a likely contribution 
of molecular gas to the mass density, and the inherent uncertainty in the measured pa- 
rameters used as input to our analytical model, the main conclusion from this work is 
our confirmation of th e flarin g of H I in dwarf galaxies. We also confirm the prediction by 
Brinks. Walter, fc OttI ( 12002| ) that the scale height in dwarfs is larger in absolute measure 



than in spiral galaxies. In dirrs, we find scale heights from 200-400 pc in the inner regions 
to 600-1000 pc out to the last measured point. 

Measuring the H I scale height in dwarf galaxies directly is fraught with difficulties, even 
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in the relatively ideal case of an edge-on dirr. The projected Hi distribution on the sky 
is a combination of a flaring gas layer seen edge-on and a possible warp. It bec omes a 



non-trivia! 



1995 



1996 



task to extrac t the intrinsic gas density and scale height distribution (lOUing 



Warmels 



1988 



Roychowdhury et al. 



20101 ). In the case of more face-on dwarfs, 



a break in the spectral correlation functio i i can be used to in: 



( lElmegreen. Kim, fc Staveley-Smith 



2001 



Padoan et al. 



'er the scale height of the gas 



20011 ) but although this has been 



shown to work in principle in the LMC, this method requires high spatial resolution, high 
sensitivity data which is beyond what can routinely be obtained with current instruments. 

As an independent check on our results, we compared the derived radial distribution of 
the H I scale height with the size distribution of type 3 H I holes, i.e., roughly spherical cavities 
in the Hi layer created by the combined effec ts of stellar winds and supernovae in young 



stellar associations (see 



Bagetakos et al. 



2010 



for details). These structures are expanding 



and eventually their radius can exceed the thickness of the layer they find themselves in, at 
which moment they would be classified as a type 2 or type 1 hole, a type 1 hole corresponding 
to a complete blow-out. The maximum size of the type 3 holes encountered is therefore a 
lower limit on the local thickness of the gas layer. A comparison of the size distribution of 
H I holes in the four sample galaxies reveals that of the 20 type 3 holes all of them have radii 
that are in agreement with them being still fully contained within the gas layer. We note 
that if we were to adopt for Ho II and IC 2574 the models with constant velocity dispersion, 
2 out of 8 holes would exceed the thickness of the gas layer in Ho II and 1 out of 5 in IC 2574, 
lending confidence to our approach and choice of input parameters. 

In addition to the uncertainties in the input parameters and how they propagate, we 
made a couple of assumptions that should be mentioned here. In our approach, the stellar 
scale height is assumed to be constant in order to obtain the vertical dispersion (Sect. H]); 
the latter is used as input parameter in our model. Conversely, if the observed values for the 
stellar dispersion were known, we could obtain the resulting stellar scale height distribution 



ah initio. For spiral galaxies , the stellar sea 



3 within the optica! 



de Griis et all (119971 ). 



radius (INarayan &: Jog 



e heig ht is shown to fiare by a factor of 2- 



2002a|) as was also shown observationally by 



Also, the error incurred in replacing {vrotjuz each z by {vrot)obs-, which is higher in 



general, as per the approximations made in Sect. 12.21 results in a slight overestimation of 
the scale height. However, the error can not be more than a few percent as the effect of 
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including the radial term itself is found to be less than ten percent and obs ervations have 
shown that {vrot)obs R,z varies by ten percent or so along z ( iKamphuisI Il993l ) . 

In our analysis we assumed a pseudo-isothermal halo. In principle, if the H I scale height 
can be determined reliably in an independent manner, our method could be expanded to 
include as free parameter the flattening of the dark matter halo. 

Finally, in Figure |5] we show a composite plot for the radial distribution of the mid- 
plane gas volume density for the four galaxies based on the exponential fits to the relations 
in Fig. [H For reference, a volume density of O.OIMqPc"^ corresponds to 0.4 H-atomcm~^. 
The flaring of gas results in a lower mid-plane density. It would be interesting to investigate 
if this decrease in volume density explains the low star formatio n efficiency, a topi c that 



will be dealt with as part of the LITTLE THINGS project (see 
description of LITTLE THINGS). 



Hunter et al. 



2011 



for a 



6 CONCLUSIONS 

In this paper we derive the H I scale height as a function of radius for a group of four dwarf 
galaxies. We use a two-component model of gravitationally-coupled stars and gas in the 
force field of a dark matter halo to obtain the Hi scale height. We extended our earlier 
analysis by explicitly addressing the issue of the rising rotation curve of dwarf galaxies. Our 
approach gives more realistic results than earlier models in which the self-gravity of the gas 
was ignored. This is especially relevant for dwarf galaxies where stars and gas dominate the 
disk gravity, roughly at par with each other. We find that Ho II has a thick H I disk through- 
out; the other three dwarf galaxies studied in this paper show fiaring Hi as a function of 
radius. Moreover, the overall H I scale height of dwarfs is larger than in spiral galaxies. This 
is due to a combination of a lower mass density in dwarf discs as a function of radius for the 
same velocity dispersion as encountered in larger spiral galaxies. The fiaring of H I in the 
outer galaxy (beyond R = ?> — ) is mainly controlled by the dark matter halo contribu- 
tion of the net potential of the galaxy. The derived gas scale height is in agreement with the 
size distribution of type 3 Hi holes, i.e., roughly spherical, expanding cavities resulting from 
the energy deposited by rapidly evolving young stars. This type of holes is still contained 
within the gas layer and hence represent a lower limit to the gas scale height. 
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Table 1. Input parameters of the different galaxies. 



Galaxy 


S*(0) 


Rd 


R20 


a"' 


PO 


Rc 




MqPC-2 


kpc 


kpc 


kms~^ 


M0 pc"^ 


kpc 


DDO 154 


5.7 


0.8 


1.22 


8 


0.028 


1.34 


Ho II 


27.8 


1.2 


3.26 


i) 7 


0.009 


10.0 










ii)14. - 0.91R if R < 3.3 














9.9 - 0.76(R - 3.3) otherwise 






IC 2574 


24.6 


2.1 


7.53 


i) 7 


0.004 


7.23 










ii)10.8 if R < 6.4 














10.8 - 0.95(R - 6.44) otherwise 






NGC 2366 


10.5 


0.5 


2.18 


9 


0.035 


1.36 



The models for Ho II and 



IC 2574 that were run for a constant gas velocity dispersion used = 7 km s ^ . 
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Table 2. Rmax (denoting the outermost extent of the Hi disk) and Rgd (the radius where the disk surface density falls to 4 
Mq pc~^) for the different galaxies 



Galaxy 


Rmax 


Rsd 




kpc 


kpc 


DDO 154 


6.8 


3.3 


Ho II 


6.8 


6.2 


IC 2574 


9.3 


9.3 


NGC 2366 


6.8 


5.0 



The Rmax values have been taken from Bagetakos et al. (2010). 
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Figure 1. Hi scale height (HWHM) versus R for the four galaxies analysed here, with and without taking into account the 

effect of the radial term in the Poisson equation (sec legend in each of the panels). The drawn lines arc exponential fits to the 
calculated scale heights. The top four panels show our calculations, assuming for each galaxy a constant velocity dispersion 
as a function of radius. As explained in the text, the inclusion of the radial term in the Poisson equation which takes care of 
the rising rotation curve, results in a sightly larger scale height in general. The two bottom panels show the inferred H I scale 
height for Ho II and IC2574 when using a gradually decreasing velocity dispersion instead (see text for details). 
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Figure 2. Plots showing the Hi scale height (HWHM) versus normalised radius (R/Ruiax left-hand panel; ij/ij^d right— hand 
panel) for the galaxies, taking into account the effect of the radial term in the Poisson equation and for Ho H and IC 2574 using 
the radially d ecreasing velocity dispe rsion for the gas. Normalisation by -Rmax is with respect to the outermost extent of the 
Hi disk f from iBagetakos et al.ir201(J l. whereas normalisation by R^d is the radius where the mass surface density of gas plus 
stars reaches 4 Mq pc~'^. For clarity, rather than the individual radial values, the exponential fits to the analytical solution are 
shown. 
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Figure 3. Plots of logarithm of the H I scale height versus radius, showing the scale height increases exponentially with radius 
for all the galaxies, but in some objects with a different scale length over different radial ranges (DDO 154, and to a lesser 
extent IC2574). 
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Figure 4. Plots of the disk surface density (stars plus gas) versus radius for all galaxies. Comparison with Figure |3] confirms 
that for each galaxy, disk surface density plays a major role in regulating the gas scale height distribution. 




Figure 5. A composite plot depicting the radial distribution of the mid-plane gas volume density for all galaxies. 



